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^ '. Abstract 

Resonant oscillations in warped disks are examined in order to explain high- 
frequency QPOs and horizontal-branch QPOs in low-mass X-ray binaries. Different 
from our previous work, addressed to the same subjects, we relax in this paper the as- 
sumption that the disks are isothermal in the vertical direction. That is, the pressure, 
^ _ density, and temperature are assumed to be distributed in the vertical direction with 

I a polytropic relation, and the polytropic index changes as the disk state changes. By 

^ ! this generalization and by some subsidiaries we can qualitatively explain, within the 

l> 

O _ 

lO . star QPOs and little variations in black-hole QPOs. We consider vertical resonances 
O ■ 

' of g-mode oscillations, since they are most appropriate to explain observations. 

CL|' Key words: accretion, accrection disks — quasi-periodic oscillations — resonance 

O ' — warp — X-rays; stars 

^' 

> ■ 

"-^ . 1. Introduction 

ci I In this paper we are interested in possible origins of high-frequency quasi-periodic os- 

cillations (QPOs) in neutron-star and black-hole low-mass X-ray binaries (LMXBs), and addi- 
tionally in horizontal-branch QPOs in luminous neutron-star LMXBs (Z sources). One of the 
important characteristics of kHz QPOs in neutron-star LMXBs is that they usually occur in 
a pair, and their frequencies change with time (see, e.g., van der Klis 2000, 2004 for reviews). 
The time change occurs in such a way as their frequency ratio decreases with an increase of 
frequencies. In the case of black-hole LMXBs, high-frequency QPOs usually also occur in a 
pair, but their frequencies change little with time, the ratio being kept close to 2 : 3 (see van 
der Klis 2004 for a review). In Z sources, another important characteristic of QPOs is present. 
It is a strong correlation between kHz QPOs and horizontal branch QPOs in their frequen- 
cies (Psaltis et al. 1999). That is, the frequency of the lower kHz QPO, i^l, and that of the 
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horizontal branch QPOs, i^hbo, are correlated in each object so that i^hbo ~ O.QSi^. 

Paying attention to the closeness to 2 : 3 of the frequency ratio in high-frequency QPOs 
in balck hole LMXBs, Abramowicz and Kluzniak (2001) and Kluzniak and Abramowicz (2001) 
presented a resonance model of disk oscillations to explain the QPOs. On the other hand, the 
importance of an external deformation of disks on resonance processes has been emphasized 
by Kato (2003, 2004a, 2004b), Kluzniak et al. (2004), and Lee, Abramowicz, and Kluzniak 
(2004). Kato (2003, 2004a, 2004b) especially considered models of resonant oscillations on 
warped disks. The models were extended to a case where the warp had precession (Kato 
2005). In these warp models, the temperature is assumed to always distribute isothermally 
in the vertical direction. In real situations, however, this will not be the case. The vertical 
structure of disks will change by the difference of the disk state resulting from difference of 
the mass-accretion rate. Considering this, we assume in this paper that the pressure, density, 
and temperature distributions in the vertical direction are polytropic, and that the polytropic 
index changes with time. By adopting this model, we can easily obtain the time variations of 
the resonant oscillations as a result of changes of the polytropic index. ^ The purpose of this 
paper is to demonstrate this generalization and to apply the results to QPOs. 

In warped disks there are four types of resonant oscillations due to the possible four 
combinations of two types of oscillations (g-mode and p-mode oscillations) and two types of 
resonances (vertical and horizontal resonances). Without considering all cases, we consider here 
only a case where the time variations of QPOs can be well explained. This is the case where 
disk oscillations are g- modes, and the resonances occur through vertical motions.^ 

2. Vertical Resonances of G-Mode Oscillations in Warped Disks 

Details of the non-linear resonant oscillations on warped disks are presented by Kato 
(2003, 2004a, 2004b). An overview of the model is sketched in figure 1 of Kato (2004a). 

We consider geometrically thin relativistic disks rotating with angular velocity, Q{r). 
The epicyclic frequency on the disk is denoted by ^(r). The oscillations on geometrically thin 
disks are classified into g-mode and p-mode oscillations (see, e.g., Kato et al. 1998; Kato 2001). 
In simplified disks the oscillations are further sub- classified by the set (m, n), where m =(0, 
1, 2...) is the number of nodes in the azimuthal direction, and n =(0, 1, 2...) is the number 
related to nodes in the vertical direction. That is, n represents the number of nodes that Ur 
(the radial component of velocity associated with oscillations) has in the vertical direction. It 
is noted, however, that Uz (the vertical component of velocity associated with oscillations) has 
(n — 1) nodes in the vertical direction. In the case of n = 0, Uz = and a series of g-modes start 
from n = 1. 

^ In Kato (2005), time variations of QPOs was interpreted as a result of time change of precession of warps. 
^ In Kato (2005), horizontal resonances of g-mode oscillations were considered. 
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A warp is generally a global pattern on disks with m = n = 1. The warp may have 
precession, but here we assume that it has no precession in order to consider idealized cases. 
On a disk deformed by the warp, we superpose g-mode oscillations with arbitrary m and n 
{n is taken to be 1 or 2, later). A g-mode oscillation with the frequency uj and (m, n) has a 
relatively large amplitude, global pattern only around the radius where 

{uj - mVLf - = (1) 

is satisfied. This can be understood if the dispersion relation for local perturbations is con- 
sidered. That is, the dispersion relation shows that the region of {uj — mVlY — > is an 
evanescent region of the g-mode oscillations. In the region where (u; — mf2)^ is smaller than k^, 
on the other hand, the oscillations have very short wavelengths in the radial direction when 
geometrically thin disks are considered. 

A non-linear interaction of a g-mode oscillation with the warp produces an oscillation 
with UJ, rh, and n, where m = m ± 1 and n = n±l (these oscillations are called hereafter 
intermediate oscillations). These intermediate oscillations resonantly interact with the disk 
at the radius where the dispersion relation for these intermediate oscillations is satisfied [see 
Kato (2004b) for detailed discussions]. There are two types of resonances, corresponding two 
types of wave modes. One is resonances that occur through motions in the vertical direction 
(vertical resonances); the other is those through motions in the radial direction (horizontal 
resonances) (see, e.g., Kato 2004a,b). In this paper we are interested in vertical resonances. 
They occur around the radius where 

(u;-mfi)2-*^fi2 _Q (-2) 

is satisfied (see e.g., Kato 2004a, b, and the Appendix in this paper), where is a num- 
ber related to n and (a polytropic index specifying the vertical structure of disks, see the 
Appendix). Hereafter, we are particularly interested in cases where n = 2 (assuming n = l) and 
n = 3 (assuming n = 2). The quantities in these cases are 

^2 = 2 + ^ and *3 = 3 + ^, (3) 

as shown in the Appendix. In the limit where the disk is isothermal in the vertical direction, 
N = oo and we have \E'2 = 2 and \E'3 = 3. They are cases that have been considered in previous 
papers (Kato 2003, 2004a, 2004b). 

Combining equations (PJ) and we find that the resonances occur at radii of 

K={'^l^^-1)Q and ^=(^3/^-1)^, (4) 

for n = 2 and n = 3, respectively. After this resonance the intermediate oscillations feedback to 
the original oscillations by non-linear interaction with the original oscillations themselves. The 
results are amplification or dampening of the original oscillations, depending on the types of 
oscillations and resonances (Kato 2004b). 
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Detailed examination shows that when m = m — 1, the frequencies of the g-mode oscilla- 

1 /2 

tions that resonantly interact with the disk at k = — l)fl are mQ — k at the radius. On the 
other hand, when m = m + 1 the frequencies of the g-mode oscillations that resonantly interact 

1 /2 

with the disk at k = — 1)Q are mQ + k there (e.g., Kato 2004a). The most observable 
oscillations are those with small m's; otherwise, they are phase-mixed. Hence, the frequencies, 
u, of typical non-axisymmetric resonant oscillations are Q — k, (i.e., m = l,m = 0), Q + k, (i.e., 
m = l,m = 2), and 2Q — k, (i.e., m = 2,m = 1). As axisymmetric resonant oscillations we have 
K (i.e., m = 0,m = 1). Axisymmetric oscillations, however, will be less observable. Hence, we 
think that the above non-axisymmetric oscillations are related to the horizontal branch QPOs, 
upper and lower kHz QPOs. Considering this, we introduce 

cuh = ^ + uj^ = 2Q — K, and cuhbo = — (5) 

and examine their frequencies at the resonance radius. 

3. Frequencies and Their Variations of Resonant Oscillations 

In the limit of isothermal disks in the vertical direction, i.e., = oo, the resonance occurs 
at K = {\/2 — when h = 2 and at k = — 1)17 when n = 3. These radii are 3.62rg and 
QAGvg, respectively (e.g., Kato 2004a). The temperature distribution in the vertical direction 
in disks is, however, generally not isothermal. Furthermore, it changes with time, depending 
on the state of the disks. In figure 1, the r/rg-^! relation given by k = (\[/^/2 — \)Vt is shown for 
the range of \I' = 2.0-4.0. In the case of ri = 2, \I' (i.e., ^2) changes from 2 (for 7 = 1 or A^ = 00) 
to 2.67 (7 = 5/3 or A^ = 1.5), while \l/3 changes from 3 (7 = 1) to 5 (for 7 = 5/3) in the case of 
n = 3. The ranges of variations of '^2 and are also shown in figure 1. 

If the value of \1/ is specified, the resonant radius is obtained. Then, from equation (0), 
the frequencies of the resonant oscillations, tun, ^^l, and whbo are derived as functions of The 
relations among these frequencies are shown in figure 2 as functions of cjh- The relation between 
2i^HBO and cjh is also shown. Important points are that: i) figure 2 is free from the mass of 
the central source, if the both horizontal and vertical scales are normalized by (M/M©)"^, and 
ii) the curves in figure 2 are universal as long as g-mode oscillations are concerned. They are 
free from A^ (and \1/), and even free from the presence or absence of precession (see Kato 2005). 
That is, the frequencies of the resonant oscillations change along the curves when A^ (and thus 
^) is changed. The relation between cuh and \1/ is shown in figure 3. In general, an increase of 
^ decreases ci;h. The range of the variation of cuh by a change of \1/ is also shown in figure 2. 

4. Summary and Numerical Estimate 

The main results of this paper are summarized in figure 2. The figure should be com- 
pared with figure 2.6 of van der Klis (2004), which summarizes observational data concerning 
QPO frequencies on a frequency-frequency diagram. This comparison suggests that the reso- 
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nant oscillations specified by tun, ^l, and ujubo "well correspond to the upper and lower kHz 
QPOs and the horizontal branch QPOs in neutron-star LMXBs, respectively. Adopting these 
identification, we qualitatively explain some basic observational characteristics of the QPOs. 

The kHz QPOs in neutron-star LMXBs usually appear in a pair, and the separation 
frequency of the twin peaks decreases as the peak frequencies increase. This observational 
characteristic is derived in our model, as shown in figure 2 (see the curves of cul and ciJh)- 
In our model the frequency change of QPOs is the result of a change of the temperature 
distribution in the vertical direction. As the temperature distribution in the vertical direction 
approaches to an isothermal one [i.e., \l/2 (the case of n = 2) tends to 2.0 or \l/3 (the case of 
n = 3) tends to 3.0], the resonance radius becomes smaller (see figure 1) and the frequencies of 
resonant oscillations increase. Observations show that the frequencies of the pair QPOs increase 
with an increase of the mass-accretion rate (van der Klis et al. 1997). Hence, if our model is 
correct, it suggests that the temperature distribution in the vertical direction approaches to the 
isothermal one as the mass-accretion rate increases. Here, let us make a quantitative estimate. 
Figure 2 shows that the 3 : 2 ratio of u;h and cjl is realized when log[a;H(M/M0)] ~ 2.93, i.e., 
LJn = 851{M / Mq)~^ . This occurs for \E'3 = 3.24 (i.e., n = 3), as shown in figure 3. This gives 
= 0.08 and 7 = 1.08. The resonance radius is then at 8.33rg, as shown in figure 1. 

One of another prominent correlations among QPO frequencies in neutron-star LMXBs is 
that the frequencies of the horizontal-branch QPOs, z/rbO) are correlated with the frequencies of 
lower kHz QPOs, z/l, by i^hbo ~ 0.08z/l. In order to compare our rersults with the observations, 
the curve of the 0.08ci;L^i^H relation is drawn in figure 2. The curve crosses the curve labelled by 
whbo around \og[iUuiM / Mq)] ~ 2.46, i.e., cuh = 288(M/M0)-^ Figure 3 then shows ^3 ~ 3.66, 
which means 3/A^ ~ 0.66 or 7 ~ 1.22. This implies that the observed correlation can be explained 
if the temperature distribution in the vertical direction is, on average, around 7 ~ 1.2. From 
figure 1 we see that the resonance occurs around r = 18. Or^ when \l/3 = 3.66. As shown above, 
both cuh : i^L = 3 : 2 and cuhbo = 0.08a;L cannot be simultaneously satisfied in a rigorous sense. 
A slight larger coefficient than 0.08, say 0.09, is compatible with the ratio 3 : 2. 

In the case of black-hole LMXBs, observations show that the frequencies of pair QPOs 
change little, and their ratio is kept to be close to 3 : 2, unlike the case of neutron-star 
kHz QPOs. As discussed in the case of neutron-star QPOs, one possibility of explaining the 
observed 3 : 2 is that it represents tun : (^l- Another one, which is better, is to consider the 
resonances of n = 2 and to regard the ratio as cul : 2co'hbo- As discussed in Kato (2004b), if 
we consider the resonance at 4.0rg, uu is equal to cjl and the ratio co'L(or uu) '■ 2co'hbo is just 
3 : 2. Figure 1 shows that the resonance at A.Ovg is realized for n = 2 when \l/2 = 2.25, which 
means 1/A^ = 0.25 (i.e., 7 = 4/3) in the present model. Figure 2 (see also figure 3) shows that 
this occurs when log[a;H(M/M0)] = 3.33, i.e., un = 2.14 x 10^{M/Mq)-\ It is noted that in 
the case where the resonance occurs through intermediate oscillations of n = 2 (not n = 3), the 
frequency variation for a change of (or 7) is weak (see the variation range of \l/2 in figure 
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2). Furthermore, cjh = 2.14 x 10^{M/Mq)~^ is compatible with observed resuhs derived by 
McChntock and Remillard (2003). It is, however, unclear why the balck-hole QPOs are the 
resonances of n = 2 (n = 1), while the neutron-star QPOs are the resonances of n = 3 (n = 2). 

More quantitative comparisons of our model with observations would be premature at the 
present stage, since the frequencies of vertical resonances are sensitive to the vertical structure 
of the disks (the vertical structure is not always polytropic). In other words, if our present 
model is correct, a comparison of our results with observations would give good information 
concerning the vertical structure of disks. 

5. Discussion 

The present resonance model naturally explains some basic characteristics of the kHz 
QPOs of neutron-star LMXBs and the high-frequency QPOs of black-hole LMXBs. It is espe- 
cially noted that the observed frequency change of kHz QPOs, which is observationally related 
to a change of the mass-accretion rate, is naturally explained if the vertical structure of disks 
changes with the change of the mass-accretion rate. In this sense, the present model seems to 
be superior to a precession model of warps (Kato 2005). In the latter model, a time variation 
of precession is required to explain the time variation of the observed kHz QPO frequencies, 
and it is unclear whether such a variation of precession is generally expected theoretically. The 
latter model with precession, however, can naturally explain the observed hectohertz QPOs as 
a manifestation of the precession of warps. In the present model of the vertical resonances, 
on the other hand, there is a problem concerning excitation. The vertical resonances do not 
excite the g-mode oscillations, but rather dampen them in the limit of = oo (Kato 2004b). 
We should carefully study in the future whether our previous results concerning the stability 
of resonancs are correct and relevant, even when N y^oo. Comparisons of the characteristics of 
the present model with those of the precession model are given in table 1. 

Many QPO models have been proposed so far. Some of them connect the observed QPO 
frequencies with the characteristic frequencies of disks, such as the orbital frequency, radial and 
vertical epicyclic frequencies and others, or their combinations. For example, in their precession 
model, Stella and Vietri (1998) identify the periastron precession frequency, Q — n [see the third 
equation in with the lower frequency of the kHz QPO. An important issue in such models is 
where preferred radii to produce specific frequencies exist. Concerning this point, these models 
are classified into various types, e.g., precession model (Stella, Vietri 1998), resonance model 
(Kluzniak, Abramowicz 2001; Abramowicz, Kluzniak 2001), and beat-frequency model (Miller 
et al. 1998). Our warp model is based on hydrodynamical wave phenomena and different from 
the models mentioned above. However, if we discuss the present model in connection with 
them, a warp can be regarded as a process to select preferred radii. 

Finally, the applicability of the present model to cataclysmic variables (CVs) is briefly 
discussed. Mauche (2002) pointed out that the frequency correlation z/rbo ~ 0.08z/l can be 
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Table 1. Comparison of two resonant models of g-mode oscillations. 



Horizontal resonances (Kate 2005) Vertical resonances (present paper) 



Time variation 



change of precession change of vertical disk structure 



Hectohertz QPOs 



precession 



precession ? 



Excitation 



O (probably) 



? 



extended to CVs. This has recently been confirmed by Warner and Woudt (2004). That is, the 
DNO (dwarf novae oscillation)-QPO relation in CVs is on the line of extension of the HBO-kHz 
QPO relation in X-ray stars. Furthermore, it is known that the DNOs in CVs change their 
frequencies, accompaning harmonics (Warner, Woudt 2004). 

The disks of CVs are Keplerian and Q and k, are almost equal. Hence, at a glance, the 
resonance conditions, i.e., equations (jU, seem not to be realized anymore. This is, however, not 
the case. The observations show that DNOs usually occur in the phase of outbursts. Near the 
transition front of outbursts, the disk thickness changes sharply. In such region, the derivation 

1 /2 

of the resonance condition, k = (\l/3 — l)fl, will be inaccurate, since the assumption of slow 
radial change of H/r is involved in the derivation. If we assume, however, that the resonance 

1 /2 

condition, n = {'^^ — l)f2, is still valid, the resonances occur even in the Newtonian Kepler 
disks, if ~ 4. Since ^3 = 3 + ^3 = 4 is realized when 1/N = 1/3, i.e., 7 = 4/3. (As 
mentioned before, \E'3 = 3 for 7 = 1 and \E'3 = 5 for 7 = 5/3.) A vertical disk structure with 
= 3 will not be unrealistic. This consideration suggests that the harmonic structure of DNOs 
can be interpreted as the oscillations of 2f2 + k, Q + k, and k, at the resonant radius. Their 
frequency ratios are 3 : 2 : 1 in the Newtonian Kepler disks. The observed frequency changes 
of DNOs are interpreted as being the result of a change of the resonance radius by a change 
of the disk structure. [See also Kluzniak et al. (2005) for an interpretation of DNOs.] The 
observed correlation between DNOs and QPOs in CVs, however, cannot be explained by the 
present model. Further considerations are needed. 

Appendix 1. Vertical Oscillations of Polytropic Disks 

We employ cylindrical coordinates (r, (f, z), whose origin is at the center of the central 
object and the 2;-axis is in the vertical direction perpendicular to the disk plane. Let us consider 
a disk where the pressure, p, and density, p, are related in the vertical direction with a polytropic 
relation, i.e., p = Kp^^^^^ , where A^ is the polytropic index. Then, the vertical integration of 
the hydrostatic balance in the vertical direction gives 




(Al) 




(A2) 



and 



7 



^2 X N+l 

PG=VQQ{r)[l - — ] , (A3) 



where subscript represents the quantities in the equihbrium state, and 00 represents those on 
the equatorial plane (e.g., Kato et al. 1998). The half-thickness of the disk, H{r), is related to 
Poo and poo by 

nlH^ = 2{N + l)^. (A4) 
Poo 

We consider here a vertical oscillation of the disk, assuming, for simplicity, that it has 
no velocity in the horizontal direction. The perturbation associated with the oscillation is 
assumed to be proportional to exp[i{ujt — rmp)], where u is the frequency of the oscillation 
and m is the number of arms in the azimuthal direction. If the oscillation occurs under the 
polytropic relation given above, the vertical component of the equation of motion gives 

dhi 

i{uj — mfl)uz = (A5) 

where hi = c^pi/po and and pi are the vertical velocity and the density perturbation asso- 
ciated with the oscillation, respectively, and Cg = 7Po/po with 7 = 1 + 1/N. The equation of 
continuity is written as 

d 

i{uj -mn)pi + —(pqU;,) = 0. (A6) 

Elimination of Uz from the above two equations gives an equation for hi as 
1 d f dhi\ (cu-mfi)2 

which is reduced to 

d^hi 2Nx dhi ^ 2N {u-mn)'' ^^_^ 
dx^ 1 ~ x'^ dx 1 — 

by changing the independent variable from z to a dimensionless variable x defined hj x = z/H. 
This equation should be solved with a boundary condition that the Lagrangian change of 
pressure is zero at the disk surface to obtain eigen- values cu's. 

The solutions of equation ()A8jl are easily obtained by expressing hi in terms of a finite 
power series of x. In the fundamental mode of oscillations, hi c(.x and {u — mQ)"^ = Q"^. In the 
first overtone, /ii oc 1 — (1 + 2A^)x^ and {u — mQ)'^ = {2 + l/N)Q'^. In the second overtone we have 
hi oca; — (l + 2A^/3)a;^ and {u — mQ)^ = (3 + 3/A^)r2^. In this way we obtain the eigen-frequencies 
of vertical oscillations as 

{tu-mny = ^M'^, (A9) 
where for the fundamental [n = 1), first overtone (ra = 2), and second overtone (n = 3) we have 



^1 = 1, (AlO) 
1 

iV' 



^2 = 2 + 4, (All) 
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and 

^3 = 3 + ^. (A12) 

For 7 = 5/3 and 7 = 4/3, we have 1/A^ = 2/3 and 1/3, respectively. In the case of isothermal 
disks, 1/A^ = 0. 
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Fig. 1. The r/rg-'^ relation obtained by solving the resonance condition k = ['^^/'^ — The disk is 
Keplerian with the Schwarzschild metric. The ranges of variation of and ^'3 arc shown. 
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Fig. 2. Frequencies (wh, whbOi ^-^d Swhbo) a-s functions of cjh- For a comparison, the 0.08wl~<^h 
relation is also shown. The ranges of variation of '^2 and 'fa are shown near the horizontal scale. 
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Fig. 3. Vt-WH relation. The disk is Keplerian with the Schwarzschild metric. The ranges of variations of 
^2 and are also shown. 
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